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Abstract — WindSat  is  a  space-based  polarimetric  microwave  ra¬ 
diometer  designed  to  demonstrate  the  capability  to  measure  the 
ocean  surface  wind  vector  using  a  radiometer.  We  describe  a  non¬ 
linear  iterative  algorithm  for  simultaneous  retrieval  of  sea  surface 
temperature,  columnar  water  vapor,  columnar  cloud  liquid  water, 
and  the  ocean  surface  wind  vector  from  WindSat  measurements. 
The  algorithm  uses  a  physically  based  forward  model  function  for 
the  WindSat  brightness  temperatures.  Empirical  corrections  to  the 
physically  based  model  are  discussed.  We  present  evaluations  of 
initial  retrieval  performance  using  a  six-month  dataset  of  WindSat 
measurements  and  collocated  data  from  other  satellites  and  a  nu¬ 
merical  weather  model.  We  focus  primarily  on  the  application  to 
wind  vector  retrievals. 

Index  Terms — Microwave  radiometer,  ocean  surface  winds, 
polarimetric,  retrieval,  WindSat. 


I.  Introduction 

THE  first  space-based  fully  polarimetric  microwave  ra¬ 
diometer,  WindSat  [1],  was  launched  in  January  2003. 
WindSat’ s  primary  mission  is  to  provide  measurements  for  the 
evaluation  of  polarimetric  microwave  radiometry  in  retrieving 
the  ocean  surface  wind  vector.  WindSat  also  provides  measure¬ 
ments  for  retrieving  sea  surface  temperature  (Ts),  columnar 
atmospheric  water  vapor  ( [V ),  and  columnar  atmospheric  cloud 
liquid  water  ( L ). 

The  polarization  properties  of  an  electromagnetic  wave 
can  be  fully  characterized  by  measuring  the  modified  Stokes 
vector.  The  modified  Stokes  vector  includes  the  vertical  and 
horizontal  polarizations  and  the  third  and  fourth  Stokes  param¬ 
eters  [Tv.  Th.  ThyX-}].1  Modeling  and  aircraft  measurements 
have  shown  that  Tv  and  7),  are  even  periodic  functions  of 
<f)R,  and  T3  and  T4  are  odd  periodic  functions  of  4>r  [2]-[4] 
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where  <pR,  the  relative  wind  direction,  is  defined  as  the  compass 
wind  direction  minus  the  radiometer  look  direction.  Therefore, 
dual-polarization  radiometers,  which  measure  only  Tv  and  Th, 
do  not  provide  enough  information  to  unambiguously  retrieve 
wind  direction.  However,  a  fully  polarimetric  radiometer  such 
as  WindSat,  which  also  measures  T3  and  T4,  provides  sufficient 
information  to,  at  least  in  principle,  retrieve  the  ocean  surface 
wind  vector. 

Dual-polarization  observations  from  radiometers  such  as 
the  Special  Sensor  Microwave/Imager  (SSM/I)  [5]  and  the 
Advanced  Microwave  Scanning  Radiometer-EOS  (AMSR-E) 
[6]  have  been  used  to  retrieve  ocean  surface  wind  speed2 
(W),  Ts,V,  and  L  with  both  statistical  and  physically  based 
methods.  Most  statistical  regression  algorithms  empirically 
derive  regression  coefficients  for  the  retrieved  parameters 
using  collocated  in  situ  measurements  or  retrievals  from  other 
satellites  (see  e.g.,  [7]-[9]).  Wentz  and  Meissner  [10]  used  a 
multiple  linear  regression  algorithm  where  the  coefficients  are 
determined  using  brightness  temperatures  (Th s)  simulated  with 
a  physically  based  model  function.  Wentz  [11]  used  physically 
based  model  functions  for  the  SSM/I  Th s,  to  obtain  a  set  of  four 
equations  in  four  unknowns  which  are  solved  using  an  iterative 
procedure.  Wentz  and  Meissner  [10]  also  outlined  a  nonlinear 
iterative  retrieval  algorithm,  but  they  do  not  discuss  retrieval 
results  obtained  with  the  algorithm. 

Algorithms  for  retrieving  wind  direction  from  a  polarimetric 
microwave  radiometer  have  previously  been  investigated  using 
aircraft  measurements  and  simulated  data.  Studies  using  aircraft 
measurements  are  limited  in  scope  because  only  a  relatively 
small  amount  of  data  is  available.  Piepmeier  and  Gasiewski  [4] 
used  aircraft  data  to  retrieve  wind  direction  using  measurements 
of  T,..  T), .  andT3  at  10.7  and  37  GHz  and  Tv  and  7),  at  18.7  GHz. 
Their  algorithm  used  maximum-likelihood  estimation  (MLE) 
for  separate  retrievals  of  wind  direction  and  wind  speed  and  it¬ 
erated  between  the  two  retrievals  to  arrive  at  a  final  wind  vector 
solution.  Liu  and  Weng  [12]  used  simulated  polarimetric  data 
to  demonstrate  wind  vector  retrievals  using  a  physical  inversion 
method.  Their  retrieval  algorithm  used  the  polarimetric  mea¬ 
surements  for  only  one  frequency. 

WindSat  provides  the  first  opportunity  to  evaluate  the 
wind  vector  retrievals  from  polarimetric  radiometer  data  on  a 
global  scale.  We  previously  described  an  empirically  derived 
combined  statistical  and  MLE  algorithm  for  retrieving  ocean 
surface  wind  vectors  from  WindSat  measurements  [13].  Here 

throughout  this  paper  wind  speed  refers  to  the  equivalent  neutral-stability 
wind  speed  at  a  10-m  reference  height. 
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TABLE  I 

Nominal  NEDT  Values  for  Rain-Free  Ocean 
Retrieval  Cells  (in  Kelvin) 


Freq.,  GHz 

6.8 

10.7 

18.7 

23.8 

37 

Tv,  Th 

t3,  t4 

0.111 

0.036 

0.053 

0.033 

0.049 

0.063 

0.023 

0.036 

we  present  a  nonlinear  optimization  algorithm  which  simulta¬ 
neously  retrieves  W.  4>r ,  Ts -  V.  and  L  from  WindSat  7},s.  The 
algorithm  is  designed  to  produce  four  solutions  (ambiguities) 
for  each  set  of  WindSat  7),s  using  a  physically  based  forward 
model. 

We  begin  by  describing  in  Section  II  the  WindSat  dataset  used 
in  our  study.  In  Section  III  we  describe  our  retrieval  algorithm. 
Section  IV  describes  the  first  generation  geophysical  forward 
model  function  used  in  our  retrieval  algorithm.  We  then  present 
retrieval  performance  results  and  a  discussion  of  the  results  in 
the  final  two  sections. 

II.  WindSat  Data  Description 

A  description  of  the  WindSat  sensor  and  data  processing 
system  is  provided  in  [  1  ] .  The  data  processing  system  produces 
sensor  data  records  (SDRs)  which  contain  7),s,  geolocation  in¬ 
formation  and  data  quality  information  for  1 6  separate  channels 
at  five  different  frequencies.  Tj,s  for  the  full  modified  Stokes 
vector  are  provided  at  10.7,  18.7,  and  37  GHz.  Dual-polariza¬ 
tion  measurements,  vertical  and  horizontal,  are  provided  at  6.8 
and  23.8  GHz.  The  antenna  temperature  measurements,  which 
have  a  different  beamwidth  and  relative  pointing  angle  for  each 
frequency,  are  resampled  and  averaged  to  provide  collocated 
Tf,s  for  the  SDRs  at  a  common  resolution.  The  SDRs  used  for 
this  paper  have  an  effective  field  of  view  (EFOV)  or  footprint 
of  approximately  40  km  x  60  km.  The  nominal  effective  noise 
equivalent  differential  tempertures  (NEDT)  for  ocean  scenes 
after  resampling  and  beam  averaging  were  given  in  [13]  and 
are  repeated  here  in  Table  I. 

WindSat  was  designed  with  both  one-look  and  two-look  ca¬ 
pability  with  measurements  taken  in  both  the  forward  and  aft 
viewing  directions.  The  width  of  the  forward  swath  is  about 
950  km  and  the  width  of  the  aft  swath  is  about  350  km  where 
the  swath  width  is  defined  to  be  the  arc  length  on  the  Earth’s 
surface  where  there  are  common  measurements  available  for  all 
WindSat  frequencies  (except  6.8  GHz  due  to  the  6.8-GHz  horn 
position  on  the  edge  of  the  swath).  There  are  80  pixels  in  the  for¬ 
ward  scan  with  an  approximate  spacing  of  12.5  km  along  scan 
and  along  track.  Retrievals  are  performed  for  all  80  pixels  in  the 
forward  scan.  The  common  swath  with  6.8-GHz  measurements 
contains  63  pixels  but  due  to  rolloffs  at  the  edge  of  the  swath 
in  the  Tb s  at  6.8  GHz  we  only  use  the  6.8-GHz  measurements 
from  55  pixels  in  the  retrievals.  Due  to  the  narrowness  of  the 
aft  swath  our  initial  forward  modeling  and  retrieval  efforts  have 
focussed  on  the  forward  swath.  The  discussion  in  the  remainder 
of  this  paper  applies  only  to  the  SDRs  and  retrievals  from  the 
forward  swath. 

The  SDRs  used  for  the  retrieval  studies  described  in  this  paper 
were  produced  with  version  1.8.1  of  the  WindSat  ground  data 
processing  system.  We  use  six  months  (September  1,  2003  to 


February  28,  2004)  of  WindSat  SDRs  using  every  third  day  for 
retrieval  analysis  with  the  remainder  reserved  as  training  data  to 
develop  empirical  corrections  to  the  geophysical  forward  model 
(as  explained  in  Section  IV).  Data  are  excluded  if  the  7},s  are 
outside  of  physical  bounds  for  ocean  scenes  or  the  Earth  inci¬ 
dence  angles  are  more  than  0.5°  from  their  nominal  values.  Data 
are  also  excluded  for  rain,  ice,  radio-frequency  interference  at 
10.7  GHz  [14],  land  contamination,  inland  lakes,  for  satellite 
attitude  anomalies  and  if  less  than  60%  of  the  measurements 
nominally  used  for  beam  averaging  are  available.  Rain  is  as¬ 
sumed  to  be  present  if  the  retrieved  cloud  liquid  water  is  greater 
than  0.2  mm.  A  more  conservative  rain  flag,  which  is  based  on 
a  flag  developed  for  SSM/I  [15],  was  used  for  the  training  set. 
Rain  was  considered  to  be  present  if  any  of  the  following  con¬ 
ditions  were  satisfied: 

T37  -  0.979Th37  <  55 
1.175T„18'7  -  30  >  T37 
T^18'7  >  170 
Tf  >  210. 

III.  Retrieval  Algorithm 
A.  Optimal  Estimation 

Our  retrieval  algorithm  uses  an  optimal  estimator  [16]  which 
is  a  Gauss-Newton  iterative  method  with  a  priori  constraints. 
The  method  is  equivalent  to  minimization  of  the  cost  function 

(x  -  tafS^ix  -  Xa)  +  X2  (1) 

where  the  superscript  T  indicates  the  matrix  transpose  and  x  is 
the  state  vector  of  quantities  to  be  retrieved.  The  state  vector 
is  comprised  of  W,c/)r,Ts,V,  and  L  or  a  subset  thereof.  The 
a  priori  constraints  are  given  by  the  a  priori  state  vector,  xa,  and 
the  a  priori  error  covariance  matrix,  Sa.  The  x2  is  a  measure 
of  the  goodness  of  fit  of  the  forward  model,  evaluated  using 
the  retrieved  state  vector,  to  the  measurements  [16],  [17].  An 
estimate  of  the  x2  is  obtained  from 

X2  =  {y-F(xi))TS-\y-F{xi))  (2) 

where  the  subscript  i  denotes  the  7 1 h  iteration.  In  this  equation, 
y  is  the  measurement  vector,  with  error  covariance  matrix  Sy. 
The  measurement  vector  is  the  set  of  WindSat  T&s  used  for  the 
retrieval.  The  state  and  measurement  spaces  are  related  through 
the  forward  model,  r/j  =  F(xi).  The  forward  model  is  described 
in  Section  IV. 

The  iteration  used  can  be  written 

Xi+1  =xa+  (S71  +  Kj  S^Ki)  1 

■Kj S~1[y-yi  +  Ki ( -  xa )] .  (3) 

The  error  covariance  of  the  solution  is  approximated  by 

Sn  =  (S-l+F^S-lKny1  (4) 

where  n  is  the  final  iteration. 


BETTENHAUSEN  et  al:  NONLINEAR  OPTIMIZATION  ALGORITHM  FOR  WINDSAT 


599 


The  matrix  K, ,  called  the  weighting  function  or  kernel,  is 
the  derivative  of  the  forward  model  with  respect  to  the  state 
parameters 


(5) 


The  kernel  is  calculated  numerically  at  each  stage  of  the  itera¬ 
tion  using  a  centered  finite  difference  scheme. 

The  convergence  criteria  we  use  compares  the  change  in  the 
state  vector  at  each  iteration  to  an  estimate  of  the  retrieval  error 
covariance  [16] 


(Xi  -  x-i ,  I  )TS-  -  Xi+i)  <  N/ 4  (6) 


where  N  is  the  dimension  of  the  state  vector. 

This  retrieval  method  requires  an  accurate  forward  model  as 
a  function  of  the  retrieval  parameters,  a  priori  estimates  for  the 
retrieval  parameters,  an  error  covariance  matrix  for  the  a  priori 
and  an  estimate  of  the  measurement  error  covariances.  At  this 
time,  we  are  using  all  of  the  WindSat  7),s  for  our  retrievals  ex¬ 
cept  Tj,  at  6.8  GHz  and  T4  at  37  GHz.  1),  at  6.8  GHz  has  been 
excluded  for  this  version  due  to  a  wind  speed  bias  problem  but 
we  expect  to  include  it  in  the  future  following  improvements  to 
our  forward  model.  T4  at  37  GHz  has  been  excluded  because 
the  signal  is  too  small  (<0.2  K)  to  improve  the  retrievals. 

We  use  constant  a  priori  for  Ts,W,V,  L.  The  a  priori  for 
Ts  and  V  are  chosen  to  be  near  the  center  of  their  respective 
ranges:  Ts  =  287  K  and  V'  =  35  mm.  The  distributions  for  W 
and  L  are  highly  skewed  so  the  a  priori  constants  are  chosen 
near  the  median  for  each  distribution:  IT'  =  7  m/s  and  L  = 
0.05  mm.  The  a  priori  error  covariance  matrix,  Sa,  is  chosen  to 
be  diagonal  with  values  that  allow  the  retrievals  to  cover  the  full 
range  of  the  retrieved  parameters.  The  square  root  of  the  a  priori 
error  covariance  values  are  12  K,  6  m/s,  50  mm,  and  1  mm  for 
Ts,  W,  V,  and  L,  respectively.  To  a  large  degree,  the  a  priori 
error  covariance  values  can  be  considered  tuning  parameters. 

Initial  analysis  of  the  retrievals  from  a  single-stage  algorithm 
with  constant  a  priori  values  showed  W  and  Ts  biases  toward 
the  respective  a  priori  values  due  to  the  large  ranges  of  values 
and  reduced  forward  model  sensitivity  at  high  W  and  low  Ts. 
Therefore,  we  use  a  two-stage  retrieval  algorithm  where  the 
first-stage  retrieval  is  performed  to  provide  more  accurate  a 
priori  values  for  the  second  stage.  In  the  first  stage,  Ts,  W,  V, 
and  L  are  retrieved  using  only  the  WindSat  7},s  for  the  vertical 
and  horizontal  polarizations.  The  forward  model  for  the  first 
stage  is  not  a  function  of  the  wind  direction.  Both  retrieval 
stages  employ  the  optimal  estimation  method  descibed  above. 
A  diagram  of  the  retrieval  process  is  shown  in  Fig.  1.  This 
two-stage  process  could  be  used  in  the  future  with  a  priori  for 
the  first  stage  derived  from  climatologies  or  numerical  weather 
models  without  biasing  our  final  retrievals.  The  primary  ad- 
vangtage  would  likely  be  faster  convergence  of  the  retrieval 
states. 

The  second  stage  retrieval  solves  for  all  five  retrieval  param¬ 
eters  and  includes  the  azimuthal  wind  direction  dependence  in 
the  forward  model.  Simultaneous  retrieval  allows  the  algorithm 


Stage  1 


Stage  2 


Final  Retrieval 


Fig.  1 .  Simplified  flowchart  of  the  retrieval  algorithm. 


to  adjust  all  five  parameters  to  optimally  match  the  forward 
model  to  the  measured  Tj, s.  WindSat  Tj, s  for  all  Stokes  parame¬ 
ters  are  used.  The  first  stage  retrievals  for  Ts,  W.  V,  L  are  used 
as  a  priori  values  for  the  second  stage.  The  square  root  of  the  a 
priori  error  covariance  values  are  6  K,  4  m/s,  5  mm,  0.5  mm,  and 
45°  for  Ts,  W,  V.  L.  and  (j>R,  respectively.  These  values  have 
been  chosen  to  be  roughly  two  to  three  times  larger  than  the  ex¬ 
pected  root  mean  squre  (RMS)  errors  in  the  first  stage  retrievals 
(the  a  priori  values).  As  with  the  first  stage  a  priori  error  co- 
variance,  these  values  can  be  treated  as  tuning  parameters.  The 
method  used  to  obtain  the  wind  direction  a  priori  is  described 
below. 

The  procedure  outlined  by  (3)— (6)  is  equivalent  to  minimiza¬ 
tion  of  the  cost  function  (T),  as  noted  above.  The  x'2  will  have 
multiple  local  minima  primarily  due  to  the  dependence  of  the 
forward  model  on  </>r\  so,  the  choice  of  a  priori  may  effectively 
choose  a  local  minima.  Multiple  solutions  or  “ambiguities”  can 
therefore  be  obtained  by  performing  separate  retrievals  for  mul¬ 
tiple  a  priori  state  vectors.  We  use  four  a  priori  state  vectors, 
where  only  the  wind  direction  differs,  to  obtain  four  ambigu¬ 
ities.  The  first  a  priori  wind  direction  is  obtained  from  the  arc 
tangent  of  the  first  harmonic  upwind  and  crosswind  components 
of  the  wind  vector  from  the  two-stage  regression  algorithm  de¬ 
scribed  in  [13].  The  three  additional  a  priori  wind  directions  are 
chosen  to  be  90°,  180°,  and  270°  from  the  regression  result.  The 
retrieved  ambiguities  are  ranked  by  the  corresponding  x2  with 
the  first  rank  ambiguity  having  the  lowest  x2- 

The  remaining  information  needed  by  the  retrieval  algorithm 
is  an  estimate  of  the  measurement  error  covariance  matrix,  Sy. 
We  include  both  the  effects  of  measurement  noise  and  forward 
model  errors  in  Sy.  A  diagonal  measurement  error  covariance 
matrix  is  often  used  to  simplify  the  retrieval  calculations  and  the 
task  of  estimating  Sy.  However,  a  diagonal  Sy  does  not  account 
for  correlations  in  the  errors  for  different  channels.  While  the 
measurement  noise  in  different  channels  is  uncorrelated,  there 
are  significant  correlations  between  the  forward  model  errors. 
There  are  many  possible  contributions  to  these  correlations — we 
provide  two  examples  here.  The  atmospheric  absorption  and 
emission  are  largely  unpolarized  so  that  errors  in  estimating 
them  will  be  common  to  all  polarizations  at  a  given  frequency. 
Forward  model  errors  for  the  same  polarization  but  different 
frequencies  may  be  correlated,  for  example,  due  to  differences 
in  the  ocean  wave  spectrum  or  sea  surface  foam  coverage  be¬ 
tween  the  global  means  (for  a  given  Ts  and  W)  and  the  ob¬ 
served  conditions.  Therefore,  including  the  off-diagonal  terms 
in  Sy  improves  the  weighting  of  the  channels  for  the  retrievals. 
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TABLE  II 

Square  Root  of  the  Diagonal  of  Sy  Used  for  the  Wind  Speed 
Range  of  7-13  m/s.  The  Superscripts  in  the  Channel 
Column  Are  the  Frequencies  in  Gigahertz 


Channel 

sJSyAi  (K) 

716.8 

X  V 

0.60 

716.8 

1h 

0.78 

rnl0.7 

±v 

0.69 

rplO.7 

1h 

0.99 

n°-7 

0.26 

T10.7 

0.09 

rplS.7 

xv 

1.02 

rplS.7 

1h 

2.02 

T18.7 

0.28 

/S 

00 

0.12 

7123.8 

1.38 

7123.8 

1h 

2.51 

rp37 

M  V 

1.76 

n7 

3.65 

Tf 

0.25 

Tf 

0.09 

This  improves  the  accuracy  of  the  state  vector  estimates  and  the 
estimates  of  x2  which,  in  turn,  improves  ambiguity  selection. 

Sy  is  estimated  using  the  differences  between  the  measured 
Tj,  s  from  the  WindS  at  SDRs  and  7/,s  simulated  with  the  forward 
model 

M 

sv  =  mJ2  (Aykf 

1  k= 1 

where 

1  M 

A Vk  =  (Vk  ~  F(xk))  -  —  ~ 

3= 1 

The  yk  are  vectors  of  the  WindSat  7),s  from  the  SDRs  and 
M  is  the  total  number  of  SDRs  used  to  estimate  Sy.  Here  the 
Xj  are  state  vectors  consisting  of  collocated  data  from  a  numer¬ 
ical  weather  model  and  other  satellite  retrievals.  The  method 
for  simulating  the  7/s  and  the  data  used  are  the  same  as  those 
used  for  empirical  corrections  to  the  sea  surface  emissivity  (Sec¬ 
tion  IV).  The  measurement  error  covariance  values  vary  signif¬ 
icantly  with  wind  speed.  For  example,  the  error  covariance  for 
the  third  Stokes  7),s  is  much  smaller  at  low  wind  speeds,  where 
the  magnitude  of  the  signal  is  small,  than  it  is  at  high  wind 
speeds.  Therefore,  to  account  for  these  variations  we  calculate 
separate  error  covariance  matrices  for  five  different  wind  speed 
ranges:  W  <  4  m/s,  4  m/s <W  <  7  m/s,  7  m/s <W  <  13  m/s, 
13  m/s  <  W  <  16  m/s,  and  II  ’  >  1 6  m/s.  The  resulting  error 
covariance  matrices  account  for  modeling  error,  measurement 
noise,  and  calibration  error.  The  estimated  covariances  also  in¬ 
clude  “matchup  noise”  due  to  spatial  and  temporal  differences 
between  the  collocated  data  used  for  simulating  7),s  and  the 
values  actually  measured  in  the  WindSat  footprint.  The  square 
roots  of  the  diagonal  elements  of  Sy  are  shown  in  Table  II  where 
the  superscripts  in  the  channel  column  are  the  frequencies  in  gi¬ 
gahertz.  Table  III  shows  the  lower  triangular  part  of  Sy  for  the 
6.8-  and  10.7-GHz  channels.  The  values  for  7),  at  6.8  GHz  are 


TABLE  III 

Sy  Terms  for  the  6.8-  and  10.7-GHz  Channels  Used  for  the  Wind 
Speed  Range  of  7-13  m/s.  Sy  Is  Symmetric,  So  Only  the  Lower 
Triangular  Part  Is  Shown.  The  Superscripts  in  the  Channel 
Column  Are  the  Frequencies  in  Gigahertz 


Sv,ij  (K2) 

716.8 

■ L  V 

bo 

7110.7 
±  V 

rl°,7 

TT0.7 

ylO.7 

716.8 

J  V 

0.36 

716.8 

±h 

0.33 

0.61 

7110.7 

±  V 

0.35 

0.39 

0.47 

7110.7 

1h 

0.38 

0.69 

0.57 

0.98 

n01 

-0.01 

-0.03 

-0.01 

-0.03 

0.07 

TW.l 

0.01 

0.01 

0.01 

0.02 

-0.01 

0.01 

included  for  information  even  though  that  channel  is  not  cur¬ 
rently  in  the  retrievals.  The  values  in  Table  III  illustrate  the  cor¬ 
relations  between  the  errors  for  the  Tv  and  Th  channels  between 
frequencies  and  polarizations.  The  values  shown  in  Tables  II  and 
III  are  used  for  the  wind  speed  range  of  7-13  m/s.  The  values 
in  Table  II  decrease  up  to  50%  for  the  lowest  wind  speed  range 
and  increase  up  to  100%  for  the  highest  wind  speed  range.  The 
diagonal  elements  of  Sy  for  the  Th.  channels  are  larger  than 
those  for  the  Tv  channels  at  the  same  frequency  because  7), 
is  more  sensivite  to  changes  in  the  wind  vector  and  the  atmo¬ 
spheric  parameters. 

B.  Median  Filtering 

We  apply  a  spatial  vector  median  filter  (MF)  to  the  retrieval 
cells  to  correct  isolated  errors  in  the  ambiguity  selection  based 
on  the  x2  ranking.  The  MF  cost  function  for  a  given  retrieval 
cell  is  computed  on  a  7  x  7  cell  box  (in  scan-based  coordinates) 
centered  on  and  including  that  cell  [18].  The  median  filter  can 
be  initialized  using  the  first  rank  retrieval  from  the  optimal  esti¬ 
mation  results.  Alternatively,  the  median  filter  can  be  initialized 
with  a  “nudged”  wind  field  where  a  background  wind  field  from 
an  external  data  source  is  used  to  select  the  first  or  second  rank 
ambiguity  closest  to  the  background  field.  Further  detail  on  the 
the  median  filter  is  given  in  [13].  We  use  spatially  interpolated 
wind  fields  from  the  National  Centers  for  Environmental  Pre¬ 
diction  Final  Analysis  (NCEP)  as  the  background  wind  field  for 
the  results  presented  in  this  paper. 

IV.  Forward  Model 

Our  overall  approach  to  development  of  a  forward  model  is  to 
begin  with  a  physically  based  radiative  transfer  model  and  then 
apply  empirical  corrections  to  better  match  the  measured  7),s. 
The  empirical  corrections  are  then  used  as  a  guide  for  improve¬ 
ments  to  the  radiative  transfer  model.  In  addition,  the  depen¬ 
dence  of  the  forward  model  on  empirically  derived  relationships 
is  reduced  relative  to  a  completely  empirically  derived  model, 
so  that  the  dependence  of  the  retrievals  on  the  training  data  used 
is  limited.  Here  we  present  an  overview  of  our  forward  model 
for  the  WindSat  7],s. 

We  have  developed  a  parameterized  forward  model  similar  to 
that  described  in  [10].  The  Th s  measured  by  the  satellite  are  the 
sum  of  the  upwelling  atmospheric  radiation,  the  reflected  down- 
welling  atmospheric  and  cosmic  background  radiation,  and  the 
direct  emission  of  the  sea  surface.  The  reflected  downwelling 
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radiation  and  the  direct  emission  are  attenuated  by  the  atmos¬ 
phere.  The  Tb s  at  each  WindSat  frequency  can  be  expressed  as 

Tv,  h  =  Tup  +  r[epTs  +  rp(ClTd  own  +  tTc )]  (7) 

^3,4  =  T~ep[Ts  -  (Tdown  +  tTc)]  (8) 

where  Tv, h  refers  to  the  7],s  for  the  vertically  or  horizontally  po¬ 
larized  Tb s,  and  T3.4  refers  to  the  7},s  for  the  third  and  fourth 
Stokes  parameters.  The  sea  surface  emissivity  for  polarization  p 
isep,  and  the  corresponding  reflectivity  is  rp  =  l—ep.  Tup  is  the 
upwelling  atmospheric  brightness  temperature  at  the  top  of  the 
atmosphere,  TaOWn  is  the  downwelling  atmospheric  brightness 
temperature  at  the  surface,  and  r  is  the  atmospheric  transmis¬ 
sivity.  Tc  is  the  cosmic  background  radiation  temperature,  which 
is  approximately  2.7  K.  The  Cl  term  is  a  correction  factor  to  ac¬ 
count  for  nonspecular  reflection  of  the  atmospheric  downwelling 
radiation  from  the  rough  sea  surface  [19].  We  are  currently  ne¬ 
glecting  the  azimuthal  dependence  of  the  reflected  downwelling 
radiation  because  it  is  difficult  to  model  empirically,  and  theoret¬ 
ical  work  in  this  area  has  been  done  only  recently  [19].  We  have 
also  neglected  the  effect  of  nonspecular  reflection  of  the  cosmic 
background  radiation  since  the  effect  is  small. 

A.  Atmospheric  Parameterization 

We  use  a  one-layer  isotropic  atmosphere  approximation  be¬ 
cause  the  WindSat  frequency  band  set  does  not  provide  the  in¬ 
formation  necessary  to  estimate  atmospheric  profiles.  In  addi¬ 
tion,  the  one-layer  atmosphere  approximation  facilitates  rapid 
evaluation  of  the  forward  model  for  the  retrieval  algorithm.  The 
atmospheric  transmissivity  is  taken  to  be 

t  =  exp[—  sec  9(A0  +  Av  +  Al)]  (9) 

where  6  is  the  Earth  incidence  angle,  and  Aq.  Ay,  and  Al  are 
the  vertical  atmospheric  absorptions  due  to  oxygen,  water  vapor, 
and  cloud  liquid  water,  respectively. 

The  upwelling  and  downwelling  atmospheric  brightness  tem¬ 
peratures  are  parameterized  in  terms  of  effective  upwelling  and 
downwelling  atmospheric  temperatures,  Tjj  and  I  f, 

Tdown  =  Tn(  1  —  t)  (10) 

Tup  =  Tu(1-t).  (11) 

Values  for  Tu,Tn,  Aq,  Ay,  and  Al  are  computed  at  each 
frequency  so  that  the  parameterized  forward  model  matches 
a  plane-parallel  atmospheric  radiative  transfer  model.  The  ra¬ 
diative  transfer  calculation  uses  the  dry-air  (primarily  oxygen) 
and  water  vapor  absorption  models  given  in  [20].  Our  current 
model  excludes  precipitating  clouds.  For  nonprecipitating 
clouds  at  WindSat  frequencies,  scattering  from  cloud  liquid 
water  is  negligible  because  the  drop  size  is  small  relative  to 
the  radiation  wavelength.  The  cloud  liquid  water  absorption 
coefficient  is  therefore  proportional  to  the  cloud  liquid  water 
content  and  given  by  the  Rayleigh  approximation  [21],  We  use 
a  double  Debye  model  for  the  dielectric  constant  of  water  [22]. 
Atmospheric  profiles  from  NCEP  were  used  for  the  radiative 
transfer  calculations.  The  profiles  were  taken  from  the  1st  and 
15th  of  each  month  between  July  2001  and  June  2002  on  a 
1°  x  1°  longitude/latitude  grid.  These  data  were  filtered  to  only 
include  grid  points  that  are  in  the  ocean  between  —65°  and 


TABLE  IV 

Coefficients  for  the  Atmospheric  Parameters 
At  Each  Windsat  Frequency 


Coef. 

6.8  GHz 

10.7  GHz 

18.7  GHz 

23.8  GHz 

37  GHz 

t>D0 

244.5 

244.8 

247.9 

250.3 

245.8 

bin 

1.666 

1.813 

2.098 

2.017 

1.998 

bn  2 

-0.02772 

-0.02966 

-0.03736 

-0.03638 

-0.03345 

bn  3 

1.588e-4 

1 .660e-4 

2.262e-4 

2.223e-4 

1.91  le-4 

buo 

-0.1516 

-0.1557 

-0.2220 

-0.1607 

-0.7691 

bin 

-1.226e-3 

-2.366e-3 

- 1 7. 14e-3 

-53.47e-3 

-21.90e-3 

boo 

1.977e-2 

2.1  l2e-2 

2.792e-2 

3.723e-2 

9.594e-2 

bo  1 

-4.230e-5 

-4.390e-5 

-5.598e-5 

-7.595e-5 

-1.958e-4 

by  0 

6.3e-5 

1.8e-4 

1.695e-3 

5.185e-3 

1.820e-3 

bv  1 

2.306e-7 

5.224e-7 

2.340e-7 

1.789e-7 

6.900e-6 

by  2 

-1.558e-9 

-3.530e-9 

-1.581e-9 

2.892e-8 

-4.662e-8 

bLO 

4.04e-5 

9.89o-5 

2.913e-4 

4.572e-4 

1.0e-3 

+65°  latitude  and  at  least  75  km  from  land.  In  addition,  we  only 
included  points  where  L  <  0.3  mm  to  exclude  possible  rain. 

We  use  least  squares  fits  to  the  radiative  transfer  modeling 
results  to  calculate  the  atmospheric  parameters  for  each  forward 


model  evaluation.  The  form  of  the  fits  are 

Td  =  bo  0  +  bmV  +  bp  2V2  +  bozVi  (12) 

Tu  =  Td  +  buo  +  buiV  (13) 

Ao  =  boo  +  boiTn  (14) 

Ay  =  (byo  +  by{V  +  by2V2)V  (15) 

Al  =  bL 0(298.8  -  1.6V)L  (16) 


where  the  bij  are  the  coefficients  derived  from  the  least  squares 
fits,  V  and  L  are  in  millimeters,  and  To  and  Tu  are  in  Kelvin. 
These  functional  forms  are  similar  to  those  used  by  Wentz 
and  Meissner  [10]  with  the  following  differences.  Wentz  and 
Meissner  fit  Ay  to  only  second  order  in  V.  They  fit  To  to  fourth 
order  in  V  and  include  an  additional  term  which  includes  a  Ts 
dependence.  We  also  have  chosen  to  use  the  V  dependence  in 
the  fit  for  Al  as  a  proxy  for  cloud  temperature,  while  Wentz  and 
Meissner  used  Tg-  We  have  intentionally  avoided  connecting 
the  atmospheric  parameterization  to  Tg  to  limit  coupling  of  the 
atmospheric  and  surface  parameters  during  our  initial  forward 
model  development.  A  separate  set  of  coefficients  is  needed 
for  each  WindSat  frequency.  The  full  set  of  coefficients  used  is 
given  in  Table  IV. 

B.  Sea  Surface  Emissivity 

It  follows  from  reflection  symmetry  properties  of  polarimetric 
scattering  and  emission  from  the  sea  surface  that  Tv  and  7),  are 
even  periodic  functions  of  fo,  and  T3  and  T4  are  odd  periodic 
functions  of  fit  [2],  A  Fourier  cosine  series  for  Tv  and  7),  and 
Fourier  sine  series  for  T3  and  T4,  expanded  to  the  second  har¬ 
monic  in  <j>R,  can  be  used  to  accurately  represent  the  sea  surface 
emission  [3],  [23],  [24], 

We  consider  the  emissivity  and  reflectivity  of  the  sea  sur¬ 
face  as  determined  by  a  combination  of  the  effects  from  large- 
scale  gravity  waves,  small  scale  capillary  waves,  and  sea  surface 
foam.  The  two-scale  model  approximation  of  dividing  the  wave 
spectrum  into  large-scale  gravity  waves  and  small-scale  capil¬ 
lary  waves  has  been  shown  to  provide  general  agreement  with 
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radiometer  brightness  temperature  measurements  from  aircraft 

[24] ,  [25],  We  use  a  two-scale  model  implementation  [19]  to  ob¬ 
tain  an  initial  approximation  of  the  sea  surface  emissivity  and 
reflectivity  and  then  use  empirically  derived  corrections  to  ac¬ 
count  for  modeling  errors  and  sea  surface  foam.  Emission  from 
sea  surface  foam  increases  with  increasing  wind  speed  due  to 
increasing  foam  coverage  of  the  sea  surface  [26],  The  presence 
of  foam  increases  the  measured  l),s  because  foam  has  a  higher 
emissivity  than  water  [27],  [28],  We  have  chosen  to  account  for 
sea  surface  foam  emission  using  only  empirical  estimates  be¬ 
cause  there  is  large  uncertainty  in  current  foam  coverage  models 
[29],  However,  this  is  a  subject  of  our  current  work  and  we  an¬ 
ticipate  that  future  versions  of  our  forward  model  will  include 
an  explicit  foam  formulation. 

We  have  used  the  Durden- Vesecky  model  of  the  sea  spectrum 

[25] ,  [30]  with  modifications  to  the  cutoff  wavenumber  between 
the  large-  and  small-scale  waves,  the  hydrodynamic  modulation 
and  ao,  the  scalar  multiplier  of  the  wave  spectra.  We  continue 
to  investigate  improvements  to  the  assumed  wave  spectrum 
model  based  on  analysis  of  WindSat  7/s  and  comparisons 
to  predictions  from  the  two-scale  model.  The  wave  spectrum 
parameters  we  are  currently  using  for  the  Durden- Vesecky 
model,  which  are  based  on  our  initial  analysis  of  WindSat  7),s, 
are  described  here.  The  scalar  multiplier,  ao,  affects  the  mag¬ 
nitude  of  all  the  harmonics  including  the  isotropic  components 
of  Tv  and  7),  where  the  magnitudes  increase  with  increasing 
ao-  We  have  used  ao  =  0.0049  as  opposed  to  the  ao  =  0.004 
used  in  [30]  and  ao  =  0.008  used  in  [24]  and  [25].  The  cutoff 
wavenumber  primarily  affects  the  magnitude  of  the  second 
harmonic  terms  where  the  magnitudes  increase  with  increasing 
cutoff  wavenumber.  We  have  chosen  a  cutoff  wavenumber 
between  the  large-  and  small-scale  waves  of  ko/2.2  versus 
ko/2.0  used  in  [30],  fco/3.375  used  in  [25],  and  fco/5.0  used  in 
[24].  The  ko  term  is  the  electromagnetic  wavenumber  for  the 
individual  WindSat  frequencies.  The  hydrodynamic  modulation 
term  is  defined  in  terms  of  the  upwind  slopes  as  in  [25]  but 
modified  to  take  on  minimum  and  maximum  values  of  0  and 
2,  respectively,  as  opposed  to  the  0.5-1. 5  range  of  [25].  This 
range  of  values  for  the  hydrodynamic  modulation  increases  the 
magnitude  of  the  first  harmonic  of  the  azimuthal  wind  direction 
dependence  of  the  third  Stokes  7],s,  which  is  indicated  by 
the  WindSat  measurements.  Our  values  for  ao  and  the  cutoff 
wavenumber  are  close  to  those  given  in  [30]  and  within  the 
ranges  of  values  in  the  previous  studies. 

We  use  a  Gaussian  model  of  the  long  wave  slope  probability 
distribution  function  in  the  two-scale  model  calculations.  The 
“modified  Stogryn’’  model  of  the  sea  water  permittivity  is  used 
[22],  with  a  fixed  sea  surface  salinity  of  34  psu.  Variations  in 
salinity  have  a  small  impact  on  Ts  retrievals  resulting  from 
small  changes  in  the  7],s  at  6.8  and  10.7  GHz.  We  plan  to  use  a 
sea  water  salinity  climatology  to  improve  sea  surface  tempera¬ 
ture  retrievals  in  future  work. 

To  develop  empirical  corrections  to  the  sea  surface  emis¬ 
sivity  we  use  values  for  the  geophysical  parameters  from  NCEP, 
QuikSCAT  [31],  SSM/1,  and  TMI  collocated  to  the  training  set 
of  WindSat  SDRs  described  in  Section  II.  We  use  Ts  values 
from  NCEP  analysis  closest  in  time  and  spatially  interpolated 
to  the  location  of  the  WindSat  SDRs.  We  use  wind  speed 


and  direction  from  QuikSCAT  retrievals  within  25  km  and 
60  min  of  the  WindSat  measurement  when  available.  Before 
collocation,  the  eight  retrieval  cells  along  both  edges  of  the 
QuikSCAT  swath  were  removed,  because  they  contain  less  than 
the  optimal  four  beam  combinations,  and  have  degraded  wind 
vectors  [32].  When  a  QuikSCAT  matchup  is  not  available  we 
use  NCEP  wind  speed  and  direction  within  1  h  of  the  NCEP 
analysis  time  and  spatially  interpolated  to  the  location  of  the 
WindSat  SDRs.  Finally,  we  use  SSM/I  and  TMI  retrievals  that 
are  averaged  into  0.25°  X  0.25°  longitude-latitude  cells  for 
V  and  L  (see  geophysical  data  at  http://www.remss.com).  The 
SSM/1  and  TMI  observations  are  collocated  to  within  25  km 
and  40  min  of  the  WindSat  observations. 

We  use  the  values  for  Ts,W,  <j)R,V,  and  L  from  the  matchup 
dataset  and  the  two-scale  model  to  calculate  values  for  emis¬ 
sivity,  e modeb  atmospheric  parameters,  and  Cl .  These  values  for 
Tup,  Tdowm  t,  and  Cl  along  with  the  Tj, s  from  the  corresponding 
SDRs  and  are  then  used  in  (7)  and  (8)  to  solve  for  the  “measured 
emissivity”  of  the  ocean  surface,  emeasured-  We  then  calculate 
empirical  corrections  to  the  emissivity  using  least  squares  fits 
to  the  difference  between  emeaSured  and  emodei  in  2-m/s-wide 
wind  speed  bins.  Our  analysis  showed  that  the  following  form 
worked  well  for  the  Tv  and  7),  channels 

^measured  Griodol  —  C/0  T"  Ql  H  Cj2Ts  T"  6'3 Tg  *  (17) 

We  neglect  variations  of  T3  and  T4  due  to  changes  in  W  and  Ts 
within  a  wind  speed  bin  and  calculate  empirical  fits  of  the  form 

emeasured  =  CjO  +  Cji  Sm((/>R)  +  Cj2  Sm(2 (f)R).  (18) 

A  different  set  of  coefficients  is  calculated  for  each  2-m/s  wide 
wind  speed  bin  and  each  WindSat  channel  where  j  denotes  the 
channel.  The  corrections  for  Tv  and  7),  are  added  to  emodel-  For 
the  T3  and  T4  channels  the  correction  is  applied  as  a  ratio  to 
the  result  from  the  two-scale  model,  e.g.,  (T  —  G/  l  /  Cj  ] , model 
and  c'-2  =  C/2/ C?2, model  where  the  model  subscript  denotes  the 
harmonic  term  calculated  from  the  two-scale  model  using  the 
mean  Ts  and  W  from  the  wind  speed  bin  and  the  nominal  Earth 
incidence  angle.  The  cjq  term  for  the  T3  and  T4  channels  is  used 
as  a  radiometer  calibration  offset. 

We  assume  that  variations  in  the  Earth  incidence  angles  that 
are  within  the  nominal  ranges  (about  0.6°)  have  a  negligible  ef¬ 
fect  on  the  corrections.  This  is  valid  for  the  Tv  and  7),  channels 
because  the  corrections  account  for  less  than  5%  of  the  total 
7),  except  at  very  high  wind  speeds.  The  effects  of  Earth  inci¬ 
dence  angle  and  sea  surface  temperature  variations  on  the  T3 
and  T4  channels  are  small  but  these  effects,  as  predicted  by  the 
two-scale  model,  are  included  since  the  empirical  corrections 
are  applied  as  a  scaling  factor.  The  corrections  do  not  account 
for  modeling  errors  in  the  azimuthal  harmonics  of  the  wind  di¬ 
rection  dependence  for  the  Tv  and  7),  channels;  therefore,  the 
wind  direction  dependence  in  the  forward  model  for  Tv  and  7), 
is  just  the  result  from  the  two-scale  model.  We  have  developed 
empirical  corrections  to  the  sea  emissivity  for  wind  speeds  up 
to  20  m/s;  therefore,  we  limit  the  results  shown  in  this  paper  to 
wind  speeds  less  than  20  m/s.  The  model  function  for  high  wind 
speeds  will  be  improved  in  future  work. 
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The  two-scale  model  is  computationally  expensive  and  the 
full  calculation  cannot  practically  be  done  during  the  retrieval 
process.  Therefore,  we  calculate  emissivity  values  with  the  two- 
scale  model  and  store  the  results  in  a  three-dimensional  lookup 
table  in  Tg,  W.  and  Earth  incidence  angle.  The  emissivity  used 
for  the  retrievals  is  calculated  by  linearly  interpolating  between 
values  in  this  table.  The  maximum  T\,  error  introduced  by  this 
interpolation  method  is  less  than  0. 1  K. 

V.  Retrieval  Performance 

Wind  vector  retrievals  are  the  primary  focus  of  this  paper. 
However,  analysis  of  the  Tg ,  V,  and  L  retrievals  provides  addi¬ 
tional  information  about  how  well  the  retrieval  algorithm  is  per¬ 
forming.  The  Tg,  V,  and  L  retrieval  performance  is  also  relevant 
to  the  wind  vector  retrievals  because  our  algorithm  simultane¬ 
ously  retrieves  all  five  geophysical  parameters.  In  this  section, 
we  present  measures  of  retrieval  performance  by  comparing  our 
WindSat  retrievals  with  NCEP  for  Tg,  QuikSCAT  retrievals  for 
wind  speed  and  direction,  and  SSM/I  retrievals  for  V  and  L 
using  the  WindSat  dataset  described  in  Section  II.  The  results 
shown  here  are  intended  to  demonstrate  the  efficacy  of  the  re¬ 
trieval  algorithm  rather  than  to  verify  accurate  calibration. 

Ambiguity  selection  has  only  a  small  impact  on  our  sea  sur¬ 
face  temperature  retrievals,  and  the  effect  on  the  water  vapor 
and  cloud  liquid  water  retrievals  is  negligible.  Therefore,  the  re¬ 
sults  we  present  in  this  section  for  those  retrievals  are  only  given 
for  the  ambiguity  selected  after  median  filtering  with  nudging. 
For  the  wind  vector  results  the  ambiguity  selection  method  is 
more  important,  and  we  consider  the  differences  that  are  related 
to  ambiguity  selection.  A  25-km  collocation  distance  window 
is  used  for  all  of  the  collocated  (matchup)  datasets.  This  dis¬ 
tance  is  roughly  half  the  diameter  of  the  WindSat  SDR  foot¬ 
print.  For  datasets  where  there  are  multiple  measurements  that 
satisfy  both  the  temporal  and  spatial  thresholds,  the  measure¬ 
ment  closest  to  the  WindSat  location  is  used.  The  retrieval  re¬ 
sults  presented  are  for  every  third  day  of  WindSat  measurements 
from  the  six-month  dataset. 

A.  Water  Vapor  and  Cloud  Liquid  Water 

We  compare  our  water  vapor  and  cloud  liquid  water  retrievals 
to  SSM/I  retrievals  (http://www.remss.com).  It  is  desirable  to 
use  the  smallest  collocation  time  window  that  is  feasible  due  to 
the  high  temporal  variability  of  cloud  liquid  water.  A  40-min 
window  is  the  smallest  window  that  can  be  used  while  still  al¬ 
lowing  for  a  sufficiently  large  dataset  with  global  coverage.  This 
yields  more  than  30  million  matchups  for  the  six-month  dataset. 
The  estimated  bias  and  RMS  errors  for  the  SSM/I  water  vapor 
retrieval  algorithm  are  0.6  and  1.0  mm,  respectively  [11].  While 
no  in  situ  measurements  for  cloud  liquid  water  over  the  ocean 
are  available  for  validation  Wentz  [11]  used  an  analysis  of  the 
distribution  of  retrieved  cloud  liquid  water  to  estimate  a  RMS 
retrieval  accuracy  of  0.025  mm. 

The  overall  differences  between  our  WindSat  water  vapor  re¬ 
trievals  and  the  SSM/I  retrievals  are  a  bias  difference  of  0.43  mm 
and  a  RMS  difference  of  1 .05  mm.  Fig.  2  shows  the  difference 
in  millimeters  between  the  WindSat  and  SSM/I  water  vapor  re¬ 
trievals  versus  the  SSM/I  water  vapor.  The  differences  were  cal- 


Fig.  2.  Difference  between  the  WindSat  and  SSM/I  water  vapor  retrievals 
versus  SSM/I  water  vapor. 


Fig.  3.  Histograms  of  the  WindSat  and  SSM/I  water  vapor  retrievals  for  more 
than  33  million  collocations. 


Fig.  4.  Difference  between  the  WindSat  and  SSM/I  cloud  liquid  water 
retrievals  versus  SSM/I  retrievals. 

culated  for  measurements  in  5-mm  water  vapor  bins.  The  max¬ 
imum  RMS  difference  is  about  2  mm  at  65  mm,  or  about  3%, 
and  the  maximum  bias  difference  is  slightly  larger  than  1  mm. 
Fig.  3  shows  the  corresponding  histograms  for  the  water  vapor 
retrievals.  The  histograms  agree  well  with  only  small  differ¬ 
ences  above  60-mm  water  vapor.  The  differences  we  show  here 
between  our  WindSat  retrievals  and  SSM/I  retrievals  are  also 
on  the  order  of  (or  better  than)  the  differences  between  various 
SSM/I  algorithms  noted  in  [33]  and  [34]. 

Fig.  4  shows  the  difference  in  millimeters  between  the 
WindSat  and  SSM/I  cloud  liquid  water  retrievals  versus  the 
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Fig.  5.  Histograms  of  the  WindSat  and  SSM/I  cloud  liquid  water  retrievals  for 
more  than  33  million  collocations. 


Fig.  6.  Difference  between  the  WindSat  Ts  retrievals  and  NCEP  Ts  results 
versus  NCEP  Ts. 

retrieved  SSM/I  cloud  liquid  water.  The  differences  were  cal¬ 
culated  for  measurements  within  0.02-mm  cloud  liquid  water 
bins.  Fig.  5  shows  the  corresponding  histograms  for  the  cloud 
liquid  water  retrievals.  The  high  spatial  and  temporal  variability 
of  cloud  liquid  water  makes  quantitative  evaluation  of  these 
results  difficult  since  a  substantial  portion  of  the  differences 
may  be  due  to  collocation  differences.  However,  the  WindSat 
and  SSM/I  retrievals  are  in  good  qualitative  agreement. 

B.  Sea  Surface  Temperature 

We  have  used  NCEP  sea  surface  temperature  data  with  an 
analysis  time  within  one  hour  of  the  WindSat  retrievals  to 
evaluate  our  sea  surface  temperature  retrievals.  The  one  hour 
window  is  small  enough  such  that  temporal  changes  in  Ts 
should  not  affect  our  results.  The  NCEP  Ts  results  are  linearly 
interpolated  in  space  to  the  WindSat  measurement  location. 
The  dataset  includes  about  25  million  matchups.  We  expect  the 
accuracy  of  the  NCEP  Ts  data  used  here  to  be  similar  to  the 
accuracy  of  the  Reynolds  optimum  interpolated  (01)  Ts  data 
[35]  which  is  produced  using  similar  methods.  The  Reynolds 
OI  Ts  data  have  an  overall  standard  deviation  error  of  about 
0.5  K  and  bias  errors  less  than  or  about  0.1  K  [36],  [37]. 

The  overall  bias  difference  between  the  WindSat  Ts  retrievals 
and  the  interpolated  NCEP  Ts  values  is  —0.12  K,  and  the  stan¬ 
dard  deviation  is  0.98  K.  Fig.  6  shows  the  difference  in  degrees 
Celsius  between  the  WindSat  retrievals  for  Ts  and  the  NCEP 
results  versus  the  NCEP  Ts  in  2°  bins.  Fig.  7  shows  the  same 


Fig.  7.  Difference  between  the  WindSat  T$  retrievals  and  NCEP  Ts  results 
versus  NCEP  wind  speed. 


Fig.  8.  Derivative  of  Tv  with  respect  to  Ts  at  6.8  and  10.7  GHz  calculated 
using  the  forward  model. 

Ts  differences  plotted  versus  NCEP  wind  speed  in  2-m/s  wind 
speed  bins.  There  is  generally  good  agreement  for  Ts  in  the  20° 
to  30  °C  range  and  for  II7  <  10  m/s.  At  lower  temperatures 
and  higher  wind  speeds  the  standard  deviation  of  the  Ts  differ¬ 
ences  increases  significantly.  There  are  two  reasons  for  the  ob¬ 
served  increases  in  the  standard  deviation  of  the  Ts  differences 
at  lower  temperatures  and  higher  wind  speeds.  One  is  primarily 
a  low-temperature  effect,  and  the  other  is  primarily  a  high  wind 
speed  effect.  It  is  difficult  to  separate  the  contribution  of  these 
effects  because  the  global  mean  Ts  decreases  with  increasing 
wind  speed. 

First,  the  Ts  retrieval  accuracy  decreases  as  Ts  decreases 
because  the  Tv  sensitivity  to  changes  in  Ts(dTv/dTg)  at  6.8 
and  10.7  GHz  decreases.  This  behavior  is  shown  in  Fig.  8 
where  dTv/dTs,  as  calculated  from  our  foward  model  for  zero 
wind  speed  and  cloud  liquid  water  and  20-mm  water  vapor,  is 
plotted  versus  Ts-  The  Ts  retrievals  primarily  rely  upon  the 
measurements  of  Tv  at  6.8  and  10.7  GHz  because  those  channels 
have  the  highest  sensivity  to  Ts  and  are  relatively  insensitive 
to  atmospheric  variations.  As  discussed  in  Section  II,  Tv  at 
6.8  GHz  is  only  used  for  55  of  the  80  retrievals  in  the  WindSat 
forward  scan.  The  7),  channels  also  have  a  significant  role  in  Ts 
retrievals  because  they  are  needed  to  separate  the  wind  speed 
and  atmospheric  contributions  to  the  Tv  measurements  from  the 
Ts  contribution.  The  variations  in  dTv/dTs  are  due  to  changes 
in  the  permittivity  of  sea  water  with  Ts  [22].  The  large  positive 
Ts  bias  near  Ts  =  0  occurs  because,  even  with  the  two-stage 
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Fig.  9.  Bias  difference  between  the  WindSat  T$  retrievals  and  NCEP  'In¬ 
versus  NCEP  wind  direction  relative  to  the  WindSat  look  direction  for  three 
wind  speed  ranges. 

retrieval  algorithm,  the  sensitivity  to  Ts  variations  is  not  strong 
enough  at  low  Ts  to  pull  the  retrieved  Ts  from  the  constant 
a  priori  value  of  287  K. 

Second,  the  Ts  retrieval  accuracy  decreases  as  wind  speed 
increases  due  to  the  increasing  dependence  of  the  Tj,s  on  the 
wind  direction.  Small  forward  model  errors  for  the  directional 
dependence  can  produce  substantial  errors  in  the  Ts  retrievals. 
Fig.  9  shows  the  mean  Ts  differences  plotted  versus  NCEP  wind 
direction  relative  to  the  WindSat  look  direction.  For  wind  speeds 
in  the  0-5-m/s  range  there  is  no  Ts  bias  with  wind  direction,  but 
there  is  a  large  bias  for  wind  speeds  in  the  range  of  10-15  m/s. 
The  difference  is  largest  near  180°  because  the  magnitude  of  the 
directional  signals  for  Tv  and  7),  are  largest  there  (since  they 
vary  as  cos  <f>R  and  cos  2 <Pr).  In  Figs.  6  and  7,  the  bias  in  Ts 
versus  wind  direction  is  averaged  over  all  wind  directions  which 
increases  the  standard  deviation  of  the  differences  but  the  mean 
is  approximately  zero. 

C.  Wind  Vectors 

We  use  QuikSCAT  wind  vector  retrievals  within  one  hour  of 
the  WindSat  measurements  to  evaluate  WindSat  wind  vector 
retrieval  performance.  This  time  window  is  chosen  as  a  compro¬ 
mise  between  minimizing  the  time  difference  and  maximizing 
global  coverage.  The  resulting  dataset  contains  more  than 
29  million  matched  retrievals.  For  comparison  to  the  results 
presented  here,  analysis  of  the  differences  between  QuikSCAT 
wind  speed  and  direction  and  in  situ  measurements  from  bouys 
[38],  [39]  show  RMS  wind  speed  differences  of  1.2  m/s.  The 
directional  differences  between  QuikSCAT  and  bouys  are  about 
20°  at  5-m/s  wind  speed  and  rapidly  decrease  to  values  in  the 
range  of  10°  to  15°  for  wind  speeds  above  about  8  m/s. 

The  histograms  of  wind  speed  and  direction  retrievals  in 
Figs.  10  and  11  show  good  agreement  between  WindSat  and 
QuikSCAT.  These  plots  are  based  on  the  selected  ambiguity 
after  median  filtering  with  nudging.  The  difference  between 
the  WindSat  and  QuikSCAT  wind  direction  histograms  near 
90°  is  due  to  differences  at  wind  speeds  below  5  m/s  where 
the  directional  signal  is  small  and,  as  a  result,  the  WindSat 
retrievals  are  less  accurate. 

Fig.  12  shows  the  difference  between  WindSat  and 
QuikSCAT  retrieved  wind  speed  versus  the  QuikSCAT  wind 
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Fig.  10.  Wind  speed  histograms  for  WindSat  and  QuikSCAT  in  1-m/s  bins  for 
more  than  29  million  collocations. 


Fig.  1 1 .  Wind  direction  histograms  for  WindSat  and  QuikSCAT  in  5°  bins  for 
more  than  29  million  collocations. 


Fig.  12.  Difference  between  WindSat  and  QuikSCAT  retrieved  wind  speeds 
in  2-m/s  wind  speed  bins. 

speeds  for  the  selected  ambiguity  after  median  filtering  with 
nudging.  The  bias  of  the  difference  is  less  than  about  0.2  m/s 
for  wind  speeds  less  than  20  m/s.  The  standard  deviation  of 
the  difference  is  below  1  m/s  for  wind  speeds  below  about 
12  m/s.  At  higher  wind  speeds  the  standard  deviation  increases 
but  remains  well  below  2  m/s.  The  overall  RMS  difference  for 
all  wind  speeds  is  0.89  m/s.  The  corresponding  overall  RMS 
differences  for  the  ambiguities  based  on  \2  rank  are  0.91,  0.97, 
1.04,  and  1.28  m/s  for  the  first  through  fourth  rank  ambiguities, 
respectively.  This  shows  that  there  are  small  but  significant 
differences  between  the  wind  speeds  retrieved  for  the  four 
ambiguities. 
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TABLE  V 

Difference  Between  WindSat  and  QuikSCAT  Wind 
Direction  Retrievals  Versus  QuikSCAT  Wind  Speed 


RMS  Wind  Direction  Difference 

W  (m/s) 

First 

MF 

MF/NG 

Closest 

2-4 

94 

82 

51 

25 

4-6 

79 

64 

37 

22 

6-8 

53 

40 

22 

15 

8-10 

31 

23 

14 

10 

10-12 

22 

17 

12 

9 

12-14 

19 

15 

12 

9 

14-16 

17 

13 

11 

9 

16-18 

16 

12 

11 

9 

18-20 

17 

12 

11 

9 

Fig.  13.  RMS  difference  between  WindSat  and  QuikSCAT  wind  direction 
retrievals  in  2-m/s  wind  speed  bins. 


The  RMS  differences  between  WindSat  and  QuikSCAT  re¬ 
trieved  wind  directions  versus  the  QuikSCAT  wind  speed  are 
shown  in  Table  V  and  Fig.  13.  The  biases,  which  are  not  shown 
in  the  table  or  the  figure,  are  all  small  (less  than  3°)  for  all  wind 
speeds  less  than  20  m/s.  In  the  table,  “First”  refers  to  the  first 
ranked  ambiguity;  “MF”  refers  to  results  with  median  filtering; 
“MF/NG”  refers  to  results  with  median  filtering  and  nudging 
with  the  NCEP  background  wind  field;  and  “Closest”  refers  to 
the  retrieved  WindSat  ambiguity  that  is  closest  to  the  QuikSCAT 
direction.  The  overall  RMS  wind  direction  difference  is  30.0° 
for  the  selected  ambiguity  (MF/NG)  and  60.2°  for  the  first  rank 
ambiguity. 

Fig.  14  shows  an  example  of  a  retrieved  wind  field  from 
WindSat  data  for  September  12,  2003  with  the  wind  speed  scale 
shown  in  the  colorbar.  For  clarity,  the  vectors  are  plotted  for  a 
subset  of  the  retrieval  cells,  but  no  averaging  is  done.  Vectors 
are  not  plotted  for  wind  speeds  less  than  3  m/s.  This  figure 
gives  qualitative  verification  that  the  wind  vector  retrieval  is 
producing  realistic  wind  fields.  The  apparent  noise  in  the  wind 
convergence  zones  is  not  unexpected  because  the  WindSat 
footprint  used  here  (40  km  x  60  km)  is  likely  too  large  to 
fully  resolve  the  fronts  and  many  of  the  abrupt  changes  in  wind 
direction  are  occurring  at  low  wind  speeds. 
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Fig.  14.  Retrieved  wind  field  from  WindSat  data  for  September  12,  2003.  The 
wind  speed  is  indicated  by  the  color.  These  retrievals  are  based  on  the  selected 
ambiguity  after  median  filtering  with  nudging  (MF/NG). 

D.  Ambiguity  Selection 

The  difference  between  the  wind  direction  retrieval  perfor¬ 
mance  for  the  first  rank  and  the  closest  ambiguities,  as  shown 
in  Fig.  13,  demonstrates  the  importance  of  ambiguity  selec¬ 
tion.  Fig.  15  shows  the  ambiguity  selection  skill  for  the  se¬ 
lected  ambiguity  after  median  filtering  with  nudging  and  for 
each  of  the  four  ranked  ambiguities.  Skill  is  defined  as  the  per¬ 
centage  of  retrievals  where  the  ranked  (selected)  ambiguity  was 
the  closest  to  the  QuikSCAT  wind  direction.  The  selected  am¬ 
biguity  (MF/NG)  is  the  closest  over  80%  of  the  time  above 
W  «  5  m/s.  The  skill  for  the  third  and  fourth  ranked  ambi¬ 
guities  is  small  above  W  =  5  m/s  so  that  the  closest  ambiguity 
is  usually  the  first  or  second  rank  ambiguity.  The  combined  skill 
for  the  first  rank  and  second  rank  ambiguities  is  always  greater 
than  50%  and  is  greater  than  90%  above  about  W  =  6  m/s. 

We  use  two  steps  for  ambiguity  selection  as  discussed  in  Sec¬ 
tion  III.  The  discussion  of  ambiguity  selection  presented  here 
focuses  on  the  first  step  of  ranking  each  ambiguity  based  on  x2- 
Detailed  evaluation  of  the  effects  of  the  second  step  of  median 
filtering,  and  in  particular  median  filtering  with  nudging,  should 
include  the  study  of  local  phenomena  such  as  at  weather  fronts 
and  the  impact  of  nudging  on  the  retrieved  wind  fields.  Such  a 
study  is  beyond  the  scope  of  this  paper. 

Fig.  15  shows  that  skill  varies  significantly  with  wind  speed. 
Skill  for  a  radiometer  also  depends  on  <j)R.  The  skill  for  the  first 
and  second  ranked  ambiguities  is  shown  in  Fig.  16  versus  the 
the  QuikSCAT  wind  direction  for  W  =  6-8  m/s.  The  first  rank 
skill  has  two  maxima  near  </;/,>  =  90°  and  270°  and  minima  near 
(pjf  =  0°  and  180°.  The  second  rank  skill  varies  inversely  with 
the  first  rank  skill.  These  features  can  be  explained  by  examining 
the  directional  dependence  of  the  7),s. 
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Fig.  15.  Ambiguity  selection  skill  versus  QuikSCAT  wind  speed. 


Fig.  1 6.  Ambiguity  selection  skill  for  (solid)  the  first  and  (dashed)  second  rank 
ambiguities  versus  the  QuikSCAT  wind  direction  for  W  =  6-8  m/s. 


Fig.  17.  Directional  dependence  of  the  T&s  (as  calculated  by  our  forward 
model)  at  8-m/s  wind  speed  and  T$  =  290  K  for  10.7  GHz. 

Fig.  17  shows  the  directional  dependence  of  the  10.7-GHz 
Tb s  (as  calculated  by  our  forward  model)  at  W  =  8  m/s  and 
Ts  =  290  K.  We  have  chosen  W  =  8  m/s  for  illustration 


purposes  because  the  wind  direction  difference  for  the  closest 
ambiguity,  as  shown  Fig.  13,  is  near  the  minumum  value  while 
there  is  a  larger  difference  between  retrieval  performance  for  the 
first  rank  and  closest  ambiguities  than  there  is  at  higher  wind 
speeds.  The  T4  signals  are  nearly  pure  second  harmonic  for  all 
frequencies  and  wind  speeds  and  are,  therefore,  approximately 
equal  (ss  0)  at  <j)R  =  0°,90°,  180°,  and  270°.  The  Ts  signals 
have  significant  first  and  second  harmonic  terms;  so,  while  Ts 
is  the  same  at  (j)R  =  0°,  180°  (T3  =  0),  T3  values  at  <j>R  =  90° 
and  270°  are  distinct.  Also,  the  ratio  of  the  first  and  second  har¬ 
monics  for  T3  varies  with  frequency  and  wind  speed,  and  there¬ 
fore,  the  <j>R  values  at  which  T->  m  0  or  has  a  local  minimum  or 
maximum  also  vary. 

The  optimal  estimation  method  effectively  weights  each  mea¬ 
sured  Tb  based  on  S~ 1  where  Sy  accounts  for  measurement 
noise  and  modeling  error.  If  we  neglect  the  off-diagonal  terms  in 
Sy,  the  measurement  and  modeling  noise  are  collectively  treated 
as  zero  mean  with  a  standard  deviation  of  y/Sy!a  (see  Table  II). 
For  W  =  8  m/s,  the  peak-to-peak  amplitude  of  the  Tv  and  7), 
directional  signals  is  less  than  the  corresponding  SVtn  values 
while  the  opposite  is  true  for  the  T3  and  T4  directional  signals. 
Therefore,  the  retrieval  of  <j>R  for  a  specified  a  priori  depends 
primarily  on  the  T3  and  T4  measurements.  However,  for  some 
<f>R,  ambiguity  selection  is  dependent  on  Tr  and  7),  .  Consider 
when  <j>R  ~  0°  or  180°.  Then  both  T3  and  T4  are  near  zero,  and 
the  T3  and  T4  contributions  to  the  \'2  will  be  the  same  for  an 
ambiguity  at  =  0°  and  a  second  ambiguity  at  <f>R  =  180°. 
Therefore,  ambiguity  selection  must  rely  on  the  Tv  and  7),  mea¬ 
surements  which  results  in  low  first  rank  skill  near  both  (f>R  =  0° 
and  180°.  Conversely,  near  4>r  =  90°  and  270°,  T4  ss  0  while 
Ts  is  negative  near  <f)R  =  90°  and  positive  near  (j>R  =  270°. 
Then  T3  and  T4  provide  sufficient  information  for  ambiguity  se¬ 
lection,  and  first  rank  skill  is  high.  The  second  rank  skill  shown 
in  Fig.  16  is  higher  than  the  first  rank  skill  near  <pn  =  40°  and 
320°,  which  indicates  a  likely  forward  model  error.  Also,  note 
that  for  the  lowest  wind  speeds  the  peak-to-peak  signal  for  the 
all  polarizations  is  less  than  the  corresponding  ^JS,ha  value, 
which  leads  to  poor  first  rank  skill  for  all  (f>R. 

We  have  chosen  to  use  four  separate  a  priori  wind  directions 
to  retrieve  four  ambiguities.  This  choice  is  supported  by  the  di¬ 
rectional  dependence  of  the  Ts  and  T4  measurements.  Fig.  18 
shows  the  retrieved  wind  direction  from  WindSat  versus  the 
QuikSCAT  wind  direction  for  all  four  of  the  ambiguities  for 
W  =  4-8  m/s  (left  panel)  and  W  =  8-12  m/s  (right  panel).  The 
contours  have  been  normalized  so  that  the  sum  of  all  points  in 
each  5°  QuikSCAT  wind  direction  bin  is  100.  This  was  done  so 
that  the  shape  of  the  overall  wind  direction  distribution,  shown 
in  Fig.  11,  would  not  obscure  the  changes  due  to  the  directional 
dependence  of  the  7/s.  There  are  four  ambiguities  for  each  <f>R 
as  can  be  seen  by  noting  four  local  maxima  along  any  vertical 
line  in  the  plots.  The  four  ambiguities  are  spaced  at  nearly  90° 
intervals  when  <pR  is  near0°,  90°,  180°,  or  270°  and  are  grouped 
in  two  pairs  when  <j>R  is  near  45°,  135°,  225°,  and  315°.  This 
pairing  of  ambiguities  near  the  local  amplitude  maxima  in  T3 
and  7~4  results  from  the  uncertainty  in  the  measurements  and 
forward  model  («  ^JSypi,  as  discussed  above).  Our  discus¬ 
sion  here  of  the  ambiguity  selection  features  has  been  simpli¬ 
fied  for  the  sake  of  illustration.  We  have  explained  the  primary 
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Fig.  18. 


Two-dimensional  histogram  of  the  ambiguities  in  5°  X 


5°  bins  of  WindSat  versus  QuikSCAT  wind  direction  relative  to  the  WindSat  look  direction. 


features  using  only  the  directional  dependence  of  the  T& s.  The 
retrieval  algorithm  is  actually  considering  the  ambiguities  in  a 
five-dimensional  space.  However,  variations  between  ambigui¬ 
ties  in  the  retrieved  values  for  the  other  parameters  are  relatively 
small. 

VI.  Conclusion 

We  have  developed  a  nonlinear  iterative  retrieval  algorithm 
for  wind  vector  retrievals  from  WindSat  data.  The  algorithm 
can  easily  be  adapted  to  use  different  subsets  of  measured 
Tb s,  and  therefore,  it  can  easily  be  adapted  for  use  with  future 
polarimetric  microwave  radiometers.  The  comparisons  of  our 
WindSat  retrievals  to  QuikSCAT  retrievals  verify  that  the 
retrieval  algorithm  is  performing  well.  The  accuracy  of  the 
retrievals  is  limited  by  measurement  noise  and  the  accuracy  of 
the  forward  model.  The  differences  between  the  forward  model 
and  the  measurements  are  currently  dominated  by  modeling 
errors  as  can  be  seen  by  comparing  the  values  given  in  Tables  I 
and  II.  The  results  presented  here  are  for  the  lowest  resolution 
footprint  of  about  40  km  x  60  km.  Work  is  ongoing  to  produce 
WindSat  7),s  and  wind  vector  retrievals  at  a  higher  resolution 
of  about  25  km  x  40  km.  For  this  higher  resolution,  beam 
averaging  will  provide  less  noise  reduction  and  measurement 
noise  will  contribute  significantly  to  the  differences  between 
the  forward  model  and  the  T3  and  T4  measurements.  It  is 
clear  that  minimizing  the  measurement  noise  in  the  third  and 
fourth  Stokes  measurements  is  important  for  measuring  wind 
direction  with  a  polarimetric  microwave  radiometer.  This  will 
be  even  more  important  as  our  forward  model  is  improved. 

The  RMS  difference  between  WindSat  and  QuikSCAT  wind 
speed  retrievals  is  less  than  1.0  m/s  for  wind  speeds  below  10 


m/s.  The  increase  in  the  standard  deviation  of  the  wind  speed 
difference  at  higher  wind  speeds  is  likely  due  to  several  factors. 
The  largest  concentration  of  WindSat-QuikSCAT  matchups 
is  in  the  midlatitudes  with  relatively  few  matchups  at  high 
latitudes  where  sustained  high  winds  are  more  prevalent.  As 
a  result,  many  of  the  high  wind  speed  cases  are  from  storms 
where  the  spatial  and  temporal  variability  of  the  wind  speeds  is 
high.  The  storm-related  localized  variabilities  would  increase 
the  wind  speed  differences  between  WindSat  and  QuikSCAT 
at  high  wind  speeds.  Another  possible  cause  of  the  increasing 
variability  at  high  wind  speeds  is  errors  in  the  forward  model 
function  for  the  directional  dependence  of  the  Tv  and  7), 
channels.  This  is  similar  to  the  effect  of  these  forward  model 
errors  on  the  Ts  retrieval  performance  that  are  shown  in  Fig.  9. 
Finally,  for  high  wind  speeds  less  training  data  are  available  for 
the  empirical  corrections  to  the  sea  surface  emissivity  which 
results  in  greater  uncertainty  in  the  forward  model. 

The  wind  direction  performance  results  are  primarily  wind 
speed  dependent.  The  magnitude  of  the  wind  direction  signal  be¬ 
comes  greater  than  the  noise  in  the  measurements  around  4  m/s. 
In  addition,  forward  modeling  errors  for  the  T3  and  T4  channels 
are  currently  on  the  order  of  0. 1  K.  The  combination  of  these  ef¬ 
fects  means  that  there  is  poor  wind  direction  performance  below 
5  m/s.  The  directional  signals  for  the  Tv  and  7),  channels  are 
superposed  on  much  larger  isotropic  signals.  The  forward  mod¬ 
eling  errors  for  the  isotropic  signals  are  on  the  order  of  0.7  K 
for  Tv  and  1.0  K  for  7),  at  10.7  GHz  and  higher  at  the  higher 
frequencies.  Therefore,  the  Tv  and  7),  directional  signals  have 
a  greater  impact  on  the  wind  direction  retrievals  above  about 
10  m/s. 

The  retrieval  error  covariance  matrix,  Sn,  calculated  from 
(4)  provides  an  estimate  of  the  expected  variance  of  the  re- 
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trievals.  The  square  root  of  the  diagonal  of  Sn  should  approxi¬ 
mate  the  standard  deviation  of  the  retrieval  errors  provided  the 
foward  model,  the  measurement  error  covariance  matrix,  and 
the  a  priori  are  good  approximations.  The  overall  mean  values 
for  the  square  root  of  the  diagonal  of  Sn  give  0.86  K  for  T5, 0.78 
m/s  for  W,  1.0  mm  for  V,  and  0.032  mm  for  L.  These  values 
are  similar  to  the  overall  standard  deviation  of  the  differences 
between  our  WindS  at  retrievals  and  collocated  data  from  NCEP, 
QuikSCAT,  and  SSM/I  as  presented  in  Section  V.  Those  values 
are  0.98  K  for  T5,  0.89  m/s  for  W,  0.95  mm  for  V,  and  0.045 
mm  for  L.  For  wind  direction,  we  computed  the  mean  values  for 
the  square  root  of  the  diagonal  of  Sn  for  the  same  2-m/s  wide 
wind  speed  bins  that  we  used  for  the  comparisons  to  QuikSCAT 
The  results  are  within  2°  of  the  values  for  differences  between 
the  closest  WindSat  ambiguity  and  the  QuikSCAT  wind  direc¬ 
tion  as  give  in  Table  V.  These  results  confirm  that  our  retrieval 
algorithm  is  performing  well. 

The  retrieval  results  described  in  this  paper  provide  a  base¬ 
line  for  WindSat  performance  and  show  that  our  retrieval  algo¬ 
rithm  effectively  retrieves  all  five  geophysical  retrieval  param¬ 
eters.  These  results  have  been  obtained  using  a  first-generation 
forward  model  function  for  the  WindSat  T& s.  We  continue  to 
improve  our  forward  model  function  and  WindSat  calibration  to 
improve  the  retrievals.  We  anticipate  that  retrieval  performance 
will  improve  significantly  as  our  forward  model  is  refined. 
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